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Proposed Running Head: The Eigenvalue Spacing 
Abstract 

The eigenvalue spacing of a uniformly chosen random finite unipotent matrix in its permu- 
tation action on lines is studied. We obtain bounds for the mean number of eigenvalues lying 
in a fixed arc of the unit circle and offer an approach toward other asymptotics. For the case of 
all unipotent matrices, the proof gives a probabilistic interpretation to identities of Macdonald 
from symmetric function theory. For the case of upper triangular matrices over a finite field, 
connections between symmetric function theory and a probabilistic growth algorithm of Borodin 
emerge. 

Key words: Eigenvalue spacing, symmetric functions, Hall-Littlewood polynomial, random matrix. 



1 Introduction 

The subject of eigenvalues of random matrices is very rich. The eigenvalue spacings of a complex 
unitary matrix chosen from Haar measure relate to the spacings between the zeros of the Riemann 
zeta function ( fll|| , 13], @])- For further recent work on random complex unitary matrices, see 



12| ], Jl6[ . The references Q and [1C] contain much of interest concerning the eigenvalues of 
a random matrix chosen from Dyson's orthogonal, unitary, and symplectic circular ensembles, for 
instance connections with the statistics of nuclear energy levels. 

Little work seems to have been done on the eigenvalue spacings of matrices chosen from finite 



groups. One recent step is Chapter 5 of Wieand's thesis [16]. She studies the eigenvalue spacings of 
a random element of the symmetric group in its permutation representation on the set {1, • • • ,n}. 
This note gives two natural ^-analogs of Wieand's work. For the first g-analog, let a £ GL(n, q) 
be a random unipotent matrix. Letting V be the vector space on which a acts, we consider the 
eigenvalues of a in the permutation representation of GL(n,q) on the lines of V. Let X s (a) be 
the number of eigenvalues of a lying in a fixed arc (l,e i27r ],0 < 8 < 1 of the unit circle. Bounds 
are obtained for the mean of X s (we believe that as n — > oo with q fixed, a normal limit theorem 
holds). A second q-analog which we analyze is the case when a is a randomly chosen unipotent 
upper triangular matrix over a finite field. A third interesting q-analog would be taking a uniformly 
chosen in GL(n,q); however this seems intractable. 

The main method of this paper is to interpret identities of symmetric function theory in a 
probabilistic setting. Section || gives background and results in this direction. This interaction 
appears fruitful, and it is shown for instance that a probabilistic algorithm of Borodin describing 
the Jordan form of a random unipotent upper triangular matrix [jjj follows from the combinatorics 
of symmetric functions. This ties in with work of the author on analogous algorithms for the 
finite classical groups ||]. The applications to the eigenvalue problems described above appear in 
Section ^. We remark that symmetric function theory plays the central role in work of Diaconis 
and Shahshahani 0] on the eigenvalues of random complex classical matrices. 



2 Symmetric functions 

To begin we describe some notation, as on pages 2-5 of §. Let A be a partition of a non-negative 
integer n = Y^i hito non-negative integral parts Ai > A2 > • • • > 0. The notation |A| = n will 
mean that A is a partition of n. Let mj(A) be the number of parts of A of size i, and let A' be 
the partition dual to A in the sense that A^ = mj(A) + mj+i(A) + • • •. Let n(A) be the quantity 
Z)i>i(* ~~ ^ i s a l so useful to define the diagram associated to A as the set of points (i, j) € Z 2 
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such that 1 < j < Aj. We use the convention that the row index i increases as one goes downward 
and the column index j increases as one goes across. So the diagram of the partition (5441) is: 



Let G\ be an abelian p-group isomorphic to @ i Cyc(pf). We write G = A if G is an abelian p-group 
isomorphic to G\. Finally, let (±) r = (1 — ~) • • • (1 - 

The rest of the paper will treat the case GL(n,p) with p prime as opposed to GL(n,q). This 
reduction is made only to make the paper more accessible at places, allowing us to use the language 
of abelian p-groups rather than modules over power series rings. From Chapter 2 of Macdonald || 
it is clear that everything works for prime powers. 

2.1 Unipotent elements of GL(n,p) 

It is well known that the unipotent conjugacy classes of GL(n,p) are parametrized by partitions A 
of n. A representative of the class A is given by 
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where Mi is the i * % matrix of the form 

/ 1 1 \ 

1 1 ••• 

1 1 ••• 



1 1 

V o o o ••• o i / 

Lemmas [l||3] recall elementary facts about unipotent elements in GL(n,p). 

Lemma 1 (j^j page 181, [Tl\]) The number of unipotent elements in GL(n,p) with conjugacy class 
type A is 

\GL(n,p)\ 

Chapter 3 of || defines Hall-Littlewood symmetric functions P\(xi,X2, •■•;*) which will be used 
extensively. There is an explicit formula for the Hall-Littlewood polynomials. Let the permutation 
w act on the x-variables by sending Xi to x w u\. There is also a coordinate-wise action of w on 
A = (Ai, • • • , A„) and 5^ is defined as the subgroup of S n stabilizing A in this action. For a partition 
A = (Ai, • • • , A n ) of length < n, two formulas for the Hall-Littlewood polynomial restricted to n 
variables are: 
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P x (x lr --,x n ;t) = [ ^ ] y w(x^---xtT\ - tXj ) 

= e -(^-^ n 

Lemma 2 T/ie probability that a unipotent element of GL(n,p) has conjugacy class of type A is 
equal to either of 



2. 



A) 

P"(i)nPA(iAA,-;i) 



Proof: The first statement follows from Lemma |] and Steinberg's theorem that GL(n,p) has 
^n(n-i) un jp t e nt elements. The second statement follows from the first and from elementary 
manipulations applied to Macdonald's principal specialization formula (page 337 of ||). Full details 
appear in ||. □ 

One consequence of Lemma Q is that in the p — > oo limit, all mass is placed on the partition 
A = (n). Thus the asymptotics in this paper will focus on the more interesting case of the fixed p, 
n — > oo limit. 

Lemma 3 

PROOF: Immediate from Lemma 0. □ 

Lemmas |] and |5| relate to the theory of Hall polynomials and Hall-Littlewood symmetric func- 
tions P|. Lemma ||, for instance, is the duality property of Hall polynomials. 

Lemma 4 (Page 181 of j^) For all partitions A ;/ u, v, 

\{G X C G\ : G\/G\ =n,G 1 = v}\ = {{d C G\ : G x /G\ =v,G x = n}\. 

Lemma 5 Let G\ denote an abelian p-group of type X, and G\ a subgroup. Then for all types [i, 
v {Jgl QG X :G 1 =HI _ 1 1 

Proof: Macdonald (page 220 of Q), using Hall-Littlewood symmetric functions, establishes for 
any partitions fi, z/, the equation: 

\{G 1 CG x :G x /G l = ii,G l = v}\ _ 1 1 
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Fixing fi, summing the left hand side over all v of size n — and applying Lemma ||] yields 
vv l{GiCG A : Gx/G! = n,Gi = v}\ _ vv |{GiCG A : Gx/ G\ = u,G 1 = ,j? s \ 

X v P^ K J lU{p)mi(X) X v P^ y ' J lU{p) mi (X) 

y^ |{GiCG A :G 1 = / i}| 
Fixing n, summing the right hand side over all v of size n — and applying Lemma |3| gives that 

E 



proving the lemma. □ 



2.2 Upper triangular matrices over a finite field 

Let T(n, p) denote the set of upper triangular elements of GL(n,p) with l's along the main diagonal. 
From the theory of wild quivers there is a provable sense in which the conjugacy classes of T(n,p) 
cannot be classified. Nevertheless, as emerges from work of Kirillov [?J and Borodin it is 
interesting to study the Jordan form of elements of T(n,p). As with the unipotent conjugacy classes 
of GL(n,p), the possible Jordan forms correspond to partitions A of n. 

Theorem [l] gives five expressions for the probability that an element of T(n,p) has Jordan form 
of type A. As is evident from the proof, most of the hard work at the heart of these formulas has 
been carried out by others. Nevertheless, at least one of these expressions is useful, and to the best 
of our knowledge none of these formulas has appeared elsewhere. P\ will denote the Hall-Littlewood 
polynomial of the previous subsection. By a standard Young tableau S of size \S\ = n is meant an 
assignment of {1, • • • , n} to the dots of the partition such that each of {1, • • • , n} appears exactly 
once, and the entries increase along the rows and columns. For instance, 

13 5 6 
2 4 7 
8 9 



is a standard Young tableau. 



Theorem 1 The probability that a uniformly chosen element ofT(n,p) has Jordan form of type A 
is equal to each of the following: 

1. p p pn (x) , where fix\(jp) is the number of complete flags of an n- dimensional 

vector space over a field of size p which are fixed by a unipotent element u of type A. 

2. ^'^n^A)' — — ^— ^— — — - , where Qh\ n (p) is a Green's polynomial as defined on page 247 of 

3. (p - l) n P x (^, j|r, jjr, • • • ; ~) E^, dimix^K^xip)' where t 1 is a partition of n, dim(x^) is the 
dimension of the irreducible representation of S n of type [A, and K^x is the Kostka-Foulkes 
polynomial. 
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(p-l) n P A (i,^-,^,-;i)cham A (p) 

4- P v p n{\) — ~ > where chain\(p) is the number of maximal length chains of sub- 
groups in an abelian p- group of type A. 

5. P\(l, \ i Jj, A, ■ • • ; -) Y^s I17=i(l — m *(A, ) ) > where the sum is over all standard Young tableaux 



o/ s/iape A, and m*(Aj) is the number of parts in the subtableau formed by {1, • • ■ ,j} which 
are equal to the column number of j. 



Proof: For the first assertion, observe that complete flags correspond to cosets GL(n,p)/B(n,p) 
where B(n,p) is the subgroup of all invertible upper triangular matrices. Note that u G GL{n,p) 
fixes the flag gB(n,p) exactly when g _1 ug € B(n,p). The unipotent elements of B(n,p) are 
precisely T(n,p). Thus the number of complete flags fixed by u is ^ p _i^ T ^ n p )\ \{d '■ 9 lu 9 € 
T(n,p)}\. It follows that the sought probability is equal to (p — l) n fix\(p) multiplied by the 
probability that an element of GL(n,p) is unipotent of type A. The first assertion then follows 
from Lemma [j]. 

The second part follows from the first part since by page 187 of ||, Q^„{p) is the number of 
complete flags of an n-dimensional vector space over a field of size p which are fixed by a unipotent 
element of type A. The third part follows from the second part and a formula for Q^ n (p) on page 
247 of ||. The fourth part follows from the third part and a formula for J2p, dim{x tl )K^\{-) in ||. 

For the fifth assertion, a result on page 197 of § gives that the number of maximal length chains 

of subgroups in an abelian p-group of type A is equal to ( ^_\l„ Ssllj=i(l m * 1 ( A ) )■ Observing 

that for a partition A of n, P\(l, |, ^, ^, ■ ■ ■ ; |) = p n P\( p -, ji, ^, ■ ■ • ; |), the result follows. □ 

As a corollary of Theorem [l], we recover the "Division Algorithm" of Borodin JIJ, which gives a 
probabilistic way of growing partitions a dot at a time such that the chance of getting A after n steps 
is equal to the chance that a uniformly chosen element of T(n,p) has Jordan type A. We include 
our proof as it uses symmetric functions, which aren't mentioned in the literature on probability 
in the upper triangular matrices. 

We remark that a wonderful application of the division algorithm was found by Borodin |lj, 
who proved asymptotic normality theorems for the lengths of the longest parts of the partition 
corresponding to a random element of T(n,p), and even found the covariance matrix. We give 



another application in Section 3.2. 



Corollary 1 (J^) Starting with the empty partition, at each step transition from a partition A to 
a partition A by adding a dot to column i chosen according to the rules 



i = 1 with probability —r, 
p 1 

i = j > 1 with probability 



p 



j v j-1 



Proof: For a standard Young tableau S, let Aj(S) be the partition formed by the entries {1, • • • ,j} 
of S. It suffices to prove that at step j the division algorithm goes from Aj_i to Aj with probability 

"5 — 3 i v i — -ml m±(A ) )? because then the probability the Borodin's algorithm gives A at 

F Aj_i(i.p.^r.^.-;p; p ^ 

step n = |A| is 

PAj (1) pi p^i ^1 " " " j p) 1 „ 1 1 1 1, 



V TT J 'p'p-'p j ' -P—u - ) = Px(i - — — ■■■■-) VTTfi - 

lip, /i 1 1 1 „m*(A,y M ' „' „2' „3> '„^Z^11^ X r,m*(A,-V 

S: S / l ape(5)=Ai=l^-iV i 'p'F'F' ' P> P PP P P S j=l P 
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as desired from part 5 of Theorem p]. The fact that the division algorithm goes from Aj-_i to 

P A .(1 

Aj with probability — — 3 , ? i p i p r — ttt(1 m*(A) ) follows, after algebraic manipulations, from 

Macdonald's principle specialization formula (page 337 of ||) 

111 1, 



P P 2 P S P i (i) r 



Kp)mi(X) 
□ 

As a remark, we observe that the division algorithm ties in with an algorithm of the author for 
growing random parititions distributed according to the n — ► oo law of the partition corresponding 
to the polynomial z — 1 in the Jordan form of a random element of GL(n,p). One version of that 
algorithm is 

Step Start with A the empty partition and N = 1. Also start with a collection of coins indexed 
by the natural numbers such that coin i has probability of heads and probability 1 — ^ of 
tails. 

Step 1 Flip coin N. 

Step 2a If coin N comes up tails, leave A unchanged, set N = N + 1 and go to Step 1. 

Step 2b If coin N comes up heads, let j be the number of the last column of A whose size was 
increased during a toss of coin N (on the first toss of coin N which comes up heads, set 
j = 0). Pick an integer S > j according to the rule that S = j + 1 with probability —rr — and 



A 

p J + 1 

S = s > j + 1 with probability ^7 rr — otherwise. Then increase the size of column S of 

P S T, S-l 



li 1 1 Mil I \ " ' ' 

1 

A by 1 and go to Step 1. 
Remarks: 

1. Probabilistic algorithms similar to the one just described for GL(n,p) were used profitably 
in to prove group theoretic results of Lusztig, Rudvalis/Shinoda, and Steinberg typically 
proved by techniques such as character theory or Moebius inversion. 

2. Observe that if one condition on the (probability zero) event that each coin comes up heads 
exactly once, the transition rules are the same as for the division lemma. 

3 Applications 

In this section we return to the problem which motivated this paper: studying the eigenvalues 
of unipotent matrices in the permutation representation on lines. Lemma || describes the cycle 
structure of the permutation action of a unipotent element a of GL(n,p) on lines in V in terms of 
the partition parametrizing the conjugacy class of a. 

Lemma 6 Let a be a unipotent element of GL(n,p) with conjugacy class of type A. Every orbit of 
the action of a on the set of lines in V has size p r for some r > 0. The number of orbits of size p r 
is 

^ P ~^h — if r>l 
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Proof:As discussed at the beginning of Section ^, the matrix a may be assumed to be 

/ M Al \ 

Mx 2 

M\ 3 ■ ■ ■ ' 

\ ••• j 

where Mi is the % * i matrix of the form 

/ 1 1 \ 

1 1 ••• 

1 1 ••• 



1 1 

V ••• 1 / 

Let Ei = Mi — Id, where Id is the identity matrix. 

From this explicit form all eigenvalues of a m ,m > are 1. Thus if a m fixes a line, it fixes it 
pointwise. Hence the number of lines fixed by a m is one less than the number of points it fixes, all 
divided by p — 1, and we are reduced to studying the action of a of non-zero vectors. It is easily 
proved that Mi has order p a , where p a_1 < i < p a . Hence if a m {x\, • • • , x n ) = (x±, • • • , x n ) with 
some X{ non-zero, and m is the smallest non-negative integer with this property, then m is a power 
of p. Thus all orbits of a on the lines of V have size p r for r > 0. 

We next claim that a pT , r > fixes a vector 

(x± , • • • , Xa x , X\ 1 -\-i , • • • , X\ 1 +\ 2 , 3?Ai+A2+l , ' ' ' , 2- A1+A2+A3 > ' ' ' ) 

if and only if 

x\ 1 +:.+\ i _ 1 + p o.+i = xa 1 +...+a,_ 1 +p"+2 = • • • = xa 1+ ...+a i; = for i : Xi > p a . 

It suffices to prove this claim when A has one part Ai of size n. Observe that the ith coordinate of 
tiP a (x\, • • • , x n ) is X)j=j {J_j) x j- If « < p a , then dP a fixes all (x\, • • • , x n ) because ( p r ) = mod p 
for r < p°. If n > p a , then a p °" fixes all (aci, • • • , x n ) such that = • • • = x n = 0, for the same 

reason. Finally, if n > p a and Xj 7^ for some j > p a , let j be the largest such subscript. Then the 
j — p a th coordinate of dP a (x\ y • • • , x n ) is equal to Xj_ p a + Xj mod p, showing that a p " does not fix 
such (xi, ■ ■ ■ ,x n ). 

This explicit description of fixed vectors (hence of fixed lines) of ctP yields the formula of the 
lemma for r > 1, because the number of lines in an orbit of size p r is the difference between the 
number of lines fixed by a pT and the number of lines fixed by a pT . The formula for the number 

of lines in an orbit of size 1 follows because there are a total of p ~} lines. □ 

p-i 

3.1 Unipotent elements of GL(n,p) 

Let a be a uniformly chosen unipotent element of GL(n,p). Each element of GL(n,p) permutes 
the lines in V and thus defines a permutation matrix, which has complex eigenvalues. Each size p r 
orbit of a on lines gives p r eigenvalues, with one at each of the p r th roots of unity. For 9 G (0, 1), 
define a random variable X 9 by letting X e (a) be the number of eigenvalues of a in the interval 
(1, e 2 ] on the unit circle. For r > 1, define random variables X r on the unipotent elements of 
GL(n,p) by 
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X r (a) = P - V - . 

p — 1 

Clearly X r (a) = if r > n. Let [y\ denote the greatest integer less than y. Lemma || implies that 

x* = Xl [0\ + x: ^' i+i+ ; ,,+Xpr lp^j • 



This relationship (analogous to one used in |16| ) will reduce the computation of the mean of X to 
similar computations for the random variables X r , which will now be carried out. 

Let E n denote the expected value with respect to the uniform distribution on the unipotent 
elements of GL{n,p). 

Theorem 2 For 1 < r < n, 

EJX r ) = — - z^lm^ ~ ^ 
p — 1 

PROOF: By Lemma |, 

^ P n (^)n p X' 1 (a)+-+X' r (a) _ p X[ («)+■■ ■+X' r _ 1 (a) 

A'jH \-\' r _ A'jH hA^-i 

Observe that 2 p r-_pr— i is the number of subgroups of G\ of type v = (r). This is because 

the total number of elements of order p r in G\ is p x '^ — p Xl ^ 1 ; and every subgroup of 

type v = (r) has p r — p r ~ 1 generators. Therefore, using Lemma ||, 

„„A P'-P" 1 v |{Gi £ G> : d = (r)}| 



'mi (A) 



/ 1 O r - » r_1 \ / 1 1 

P n (-\ P P W 



P r (l-^rr)---(l-jr) 
p — 1 

□ 

Corollary |2] uses Theorem || to bound the mean of X e . 
Corollary 2 £ n (X 9 ) = 9^ - 0(£). 

Proof: Let {y} = y — |_j/J denote the fractional part of a positive number y. Theorem ||| and the 
writing of X 9 in terms of the X r 's imply that 
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E n {x 9 ) = eE n (J2Xi)-J2{p r o}En( 



X p r-l + l + • • • + X p r 



p n 


- 1 


P 


- 1 


p n 


- 1 


P 


- 1 


p n 


- 1 


P 


- 1 



i>l r>l 



v — 1 rrt 



X p r-1 + 1 + • • • + X v 



r>l P 



r=l 



(p-l)p r (p - l)pLW P (")J- 



We suppose for simplicity that n ^ pp T + 1 for some r (the case n = pP + 1 is similar). 
Continuing, 



n -i \}og p {n)\ p r +1 +1 

P-l ^ (p- l) 2 pLiogp( ")j- r+1 (p-l) 2 n p-1 

□ 

The approach here appears to extend to the computation of higher moments, but the compu- 
tations are formidable. For example one can show that if 1 < r < s < n, then 

E n (X r X s ) = ^[^(i - ■ ■ ■ (i - p;) + E(i - ^tttiLtt) ■ • • (i - 

3.2 Upper triangular matrices over a finite field 

Let a be a uniformly chosen element of T(n,p). Recall that a is unipotent by the definition of 
T(n,p). Each element of T(n,p) permutes the lines in V and thus defines a permutation matrix, 
which has complex eigenvalues. Each size p r orbit of a on lines gives p r eigenvalues, with one at 
each of the p r th roots of unity. For 9 € (0, 1), define a random variable X 9 by letting X 9 (a) be 
the number of eigenvalues of a in the interval (1, e 27 ™ 61 ] on the unit circle. For r > 1, define random 
variables X r on the unipotent elements of T(n,p) by 

X r (a) = V - V - . 

p — 1 

Let \_y\ denote the greatest integer less than y. Lemma ^ implies that 

x° = Xl [e\ + £ x ^ 1+1+ ;" + Xpr \jro\ . 

r>l P 

As for the case of GL(n,p) this relationship will reduces the computation of the mean of X e to 
similar computations for the random variables X r . 

Let E n denote the expected value with respect to the uniform distribution on the unipotent 
elements of T(n,p). Theorem || shows that the expected value of X r is surprisingly simple. As one 
sees from the case p = 2, the result is quite different from that of Theorem [2[ 
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Theorem 3 For 1 < r < n, 

E n (X r ) = (p - 1) 

Proof: We proceed by joint induction on n and r, the base case n = r = 1 being clear. Let Prob(S) 
denote the probability that Borodin's growth algorithm yields the standard Young tableau S after 
\S\ steps. Let col(n) be the column number of n in S. With all sums being over standard Young 
tableaux, observe that 




E n { P x i+-+ x 'r) = Yl P KiS)+ - +K(S) Prob(S) 

S:\S\=n 

£ pK(s)+-+K(s) Prob{s) 

S:\S\=n,col(n)=l 

+ J2 P X '^ S)+ - +X ^ S) Prob(S) 

S:\S\=n,l<col(n)=j<r 

+ p X '^ S)+ - +K{S) Prob(S) 

S:\S\=n,col(n)>r 
£ p Al(50+.--+A;(50+lp rofe(5 / ) _L_ 



□ 



S':\S'\=n-l 



V 



+£ E /'.< s ') + '"^< s ')«p™ i -(s')(T 7 ^ T - t^t) 

i=2S':|S'|=n-l P J P 

+ E E p v l(s . )+ ... + , ;(s . )pro( , (s , )( ^ _ 

j>rS':|S"|=n-l P jV P 

p# n _i(p A 2 + -+^) +p J E„_i(/'i+- +A -i -p A 2+-+ A r) 

(p A 'i + -+^-p A 'i+- +A -i) 

(p - i) j e„_ 1 (p a 'i + - +a :- 1 ) + E n ^(p X 'l+-+ X 'r) 

(p-ir l ( n \ 



Corollary |3] follows by using Theorem || and arguing along the lines of Corollary g. 
Corollary 3 ^^"H") < £ n ( X *) < 
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